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Abstract. To investigate uncertainties in air shower simulations caused by the 
small-x regime, a model including leading-twist hard pQCD plus soft processes 
was built, which are separated by an energy dependent transverse momentum 
cut-off. We provide a fit of the cut-off to the total pp cross section for different 
PDFs using the eikonal formalism and show that for modern PDF sets there 
is only a small uncertainty in the mini-jet cross section, and hence in the final- 
state multiplicity and the number of produced muons. 
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1. Introduction 

Extensive air showers (EAS), are induced in the atmosphere due to collisons of a 
primary cosmic ray particle with air nuclei. The properties of cosmic rays, such 
as their origin, energy and composition are (partly) reconstructed from air shower 
observables. Thus, understanding the physics of the primary collision and of sub- 
sequent interactions in the atmosphere is important. The main uncertainties are 
dominated by hadronic interactions. For example, muon multiplicities at Auger, 
Kascade and CosmoLEP might hint at a poor understanding of particle pro- 
duction at ultra-high cosmic ray energies. Given that these extend far beyond those 
of terrestrial accelerators, our knowledge of hadronic interactions has to be extrap- 
olated to unknown regimes in any QCD inspired model attempting to describe the 
physics of EAS [f^. 

The QCD mini-jet cross section in pp increases rapidly with energy and domi- 
nates the total cross section already for collider energies. It is sensitive to the small-x 
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regime of perturbative QCD (pQCD), and to the break-down of the leading- twist 
approximation at small transverse momentum, which requires us to introduce a 

cut-off PTcuf 

Another source of uncertainties is low-momentum "soft" particle production 
which dominates the very forward rapidity region and hence the momentum degra- 
dation of the cosmic ray in the atmosphere. For semi-central collisions with air 
nuclei an IR-safe approach to particle production is provided by gluon saturation 
effects, see [ 3 for an application to EAS. For peripheral pA and for pp collisions 
no intrinsic semi-hard scale exists and weak-coupling methods are therefore not 
available; here, we followed instead the Dual Parton Model adjusting phe- 
nomenological parameters to (recent) accelerator data. However, our present focus 
is on uncertainties related to the perturbative QCD regime, in particular with re- 
spect to the multiplicity of secondaries (or dajet/dy) which is related to the number 
of produced muons. 

2. Mini-jet production 

The hard jet cross section in pQCD can be written as [E] 

, ,2u - E ^l/«(^l.Q^fe)^2/,(x2,Q^6)-5(^+J^^ + A:), (1) 

«2/3ay4C( PTfl b TT dt 

where the summation runs over all parton species, Xi and X2 are the light-cone mo- 
mentum fractions carried by the scattering partons, t is the momentum transfer and 
b is the impact parameter. The kinematical variables xi, X2 and t are related in a 
simple way to the rapidities y^, and px ["S*. We do not employ a phenomenolog- 
ical K-ia,ctoT since its energy dependence is not well known; this uncertainty does 
not affect our results significantly except for increasing the pt-cutoff shown below 
slightly. 

The functions f{x,Q^,b) are the PDFs. The impact parameter dependence 
is assumed to factorize, f{x,Q^,b) = T{b)f{x,Q^), with T{b) the proton overlap 
function chosen according to [ .6; . HERA has measured the PDFs over a broad range 
of = p'^, but only for 10"'*. Modern PDF packages extrapolate to lower x 
according to QCD evolution equations. We shall show results for both the (leading 
order) GRV98 and CTEQ6 parameterizations [ 7 which feature a BFKL-like growth 
of the gluon distribution at small x. To illustrate the effect of small- a; evolution, 
we also employ the old Duke-Owens (DO) parameterization, which exhibits a much 
flatter xgix, Q^) at small x. 

Due to the power-law divergence oi ^ as pt ^ ^ we need to introduce a low 
momentum cut-off. The (phenomenological) soft cross section <Tsoft = 57 mb is 
assumed to be constant at high energies. It's value is chosen to reproduce the total 
pp cross section at the lower end of relevant energies, Ecm ^ 100 GeV. 

In order to avoid partial wave unitarity violation we use the eikonal formalism [ 
IH] to define the total cross section. In this formalism the total and inelastic cross 
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sections are related to <jjet + cTsoft by 



and ainel — crtot ~ Cfel- 



3. Results 



According to S-Matrix and Regge theory the total pp cross section can be param- 
eterized in a variety of ways at high energies. A parameterization fitted to the 
available experimental data is presented in [ |U] . Using expression ^ we fitted 
PTcut to that parameterization of atot (for pp) for different PDFs. As shown in 
fig. la, for modern PDFs the cut-off increases very rapidly at high energies, which 
is related to the small-x growth of the gluon distribution. The fits for CTEQ6 
and GRV98 return very similar values for PTcuti with small discrepancies only for 
Ecm — 100 TeV. On the other hand, the DO set requires a constant cut-off within 
the relevant energy regime, which is due to the too flat small-x gluon distribution 
for this old parameterization. 

Using the above PTcut fo^ ^ given energy and PDF, we can now integrate 
over 2/4 andpT to obtain dcTjet/dy. In fig. lb we show dajet/dy for different PDFs at 
Ecm = 150 TeV, a typical primary high-energy cosmic ray. Again we observe very 
good agreement among the modern PDFs, both of which predict a much higher mul- 
tiplicity at central rapidity than the old DO set. Also, the latter gives a much wider 
rapidity plateau, implying a reduced momentum degradation in the atmosphere 
from inelastic interactions. 




(a) (b) 
Fig. 1. Fig a) Transverse momentum cut-off for leading-twist LO pQCD processes 

in pp, for various sets of PDFs. Fig b) Rapidity distribution of the (mini-) jet cross 

section in inelastic pp collisions for those PDFs. 
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4. Conclusion 

We have studied possible uncertainties in hard process extrapolated to high energies, 
which are sensitive to the small-x regime of the PDFs. We fixed pTcuti^cm) by 
fitting the total pp cross section in the eikonal formalism. While modern sets of 
PDFs with a steep increase of the gluon distribution at small x agree rather well, 
old parameterizations predict a much flatter PTcut- 

We have also shown dajet/dy for different PDFs. There is little difference 
between CTEQ6 and GRV98 both at central and forward rapidity. Since mini-jets 
dominate particle production if the soft cross section is energy independent, we do 
not expect a large uncertainty in the number of produced muons (approximately 
proportional to the multiplicity of secondary hadrons) arising from the high-energy 
part of an EAS. For DO, wich exhibts a much slower growth of the small x gluon 
density one obtains a significantly wider rapidity plateu wich would affect EAS. In 
the future, we plan to extend our model to pA collisions and implement it into an 
EAS simulation package to study directly the effects of small- a; QCD interactions 
on air shower observables. 
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